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The recently observed a-decay chains 293 - 294 j^]^-!? were produced by the fusion reac- 
tions with target '^'^^Bk and projectile ^^Ca at Dubna in Russia. The reported cross- 
sections for the mentioned reaction are a = 0.5(+l.l, — 0.4)pb and cr=1.3(+1.5,-0.6)p6 at 
E* = 35MeV and E* = 39MeV, respectively. The Q-values of a-decay and the half- 
lives LoQioTy^i^) are calculated for the a-decay chains of 293-294 j^j^-j? nuclei, within the 
framework of Isospin Cluster Model (ICM). In the ICM model the proximity energy is 
improved by using the isospin dependent radius of parent, daughter and alpha particle. 
The binding energy B{Ai, Zi) (i=l,2) of any nucleus of mass number A and atomic num- 
ber Z was obtained from a phenomenological and more genaralized BW formula given by 
[25]. The calculated results in ICM are compared with the experimental results and other 
theoretical Macro-Microscopic (M-M), RMF( with NL3 and SFU Gold forces parameter) 
model calculations. The estimated values of a-decay half-lives are in good agreement 
with the recent data. The ICM calculation is in favor of the persence of magic number at 
N = 172. 
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1 Introduction 



In the nuclear chart, the superheavy mass region are growing fast due to the availabiUty 
and advancement in the radioactive nuclear beam technology. The "^^Ca is most promi- 
nent radioactive nuclear beam at present for the synthesis of superheavy elements [!,?]• 
Flerov [3] was the first who suggested in 1969, the use of a highly neutron-rich beam of 
^^Ca for the formation of superheavy elements with neutron rich targets such as ^^"^Pu, 
'^'^^Cm and "^^"^Cf. The stabihty of these superheavy nuclei depends upon the magicity 
of the proton and neutron number either spherical and/or deformed. The spherical shell 
closure for the neutron and proton numbers are at 2, 8, 20, 28, 50, and 82 while 126 is 
only for neutron number. Since a long time, the search about the next spherical shell for 
the proton number is going on. This question, about the next doubly magic nucleus be- 
yond the '^^^Pb (Z = 82,N = 126) has attracted much attention in the nuclear structure 
physics for both the theoreticans and experimentalist. Theoretical models predict [4-8] 
that the next magicity for the proton number should occur at Z = 114, 120, orl26 and 
for the neutron number it should be at A'^ = 172 and 184. 

These superheavy nuclei undergo spontaneous decay in to successive alpha decay chains 
before spontaneous fission. The a-decay energy and half-fives of these decay chains help 
us to understand the nuclear structure of the parents as well as of daughter nuclei and 
hence gives the information about the stability of peninsula. In this paper, recently ob- 
served two isotopes of the element Z=117 with mass number 293 and 294 are studied for 
the alpha decay characteristics. These two isotopes ^^^117 and ^^^117 were produced in 
the fusion reaction between '^^Ca projectile and radioactive target '^'^^Bk nuclei [9]. The 
resulting excitation energy of the compound nucleus ^^^117 is reported E* = 39MeV and 
E* — 35MeV respectively, which favoring the 4n and 3n evaporation channel. The Isospin 
Cluster Model (ICM) calculations for the 293,294227 alpha decay chains are compared with 
the experimental results [9], the macroscopic- micrscopic(MM)[10], RMF based (using 
NL3 and SFU Gold forces parameters) [11, 12] model calculations. In our earlier work for 
the CK-decay chains calculation [13-15], preformed cluster model [16] was used and it has 
been observed that calculated results could be improved using the isospin physics in the 
model. So, the earlier used model is modified with the same idea and hence obtained very 
nice results from the present ICM calculations. 
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The 293,294]^-!^'^ alpha decay chains calculation are based on the Isospin Cluster Model 
(ICM) , which is described briefly in Section 2 and the results of calculation are presented 
in Section 3. A discussion and summary of results is given in Section 4. 

2 The Isospin Cluster Model 

The ICM model uses the dynamical collective coordinates of mass (and charge) asym- 
metry, rj = (^1 — A2)/{Ai + A2) and rjz = {Zi — Z2)/{Zi + Z2), first introduced in the 
QMFT [17,18], which are in addition to the usual coordinates of relative separation R 
and deformations (i = 1, 2) of two fragments. Then, in the standard approximation of 
decoupled R and rj motions, the decay constant A or the decay half-hfe T1/2 is defined as 

X = !^ = P,v^P. (1) 

J 1/2 

Here Pq is the cluster (and daughter) preformation probability and P the barrier pene- 
trability which refer, respectively, to the 7] and R motions. The vq is the barrier assault 
frequency. The Pq are the solutions of the stationary Schrodinger equation in 77, 

which on proper normalization are given as 

Po = ^,\^^'Kv{A,)) \'{2/A), (3) 

with i=l or 2 and i/=0,l,2,3.... Eq. (2) is solved at a fixed R = = Ct{= Ci + C2). 
Cj's taken from the Myers and Swiatecki [19] droplet model. The matter radius Cj is 
calculated as 

= c, + -^t,(^ = l,2), (4) 

where Ci denotes the half-density radii of the charge distribution and U is the neutron 
skin of the nucleus. For the U calculation Myers and Swiatecki used two parameter 
Fermi function values given in Ref.[20] and remaining cases were handled with the help 
of parameterization of charge distribution decribed below. The nuclear charge radius 
(denoted as i?oo in Ref.[21]is given by the relation: 

= 1.240>l| {1 + ^ - 0.191 X (^^^) } fm{i = 1, 2), 



where < > repersents the mean square nuclear charge radius. According to Ref.[21], 
Eq.(5) was vahd for the even-even nuclei with 8 < Z > 38 only. For nuclei with Z > 38, 
the above equation was modified by Pomorski et al.,[21] as 



OOi 



1.256A3 <^ 1 - 0.202 



/m(^ = l,2), 



(6) 



These expressions give good estimate of the measured mean square nuclear charge radius< 
>. In the present model, author have used only Eq.(5). The half-density radius, q was 
obtained from the relation: 

49 



QOi 



7 

1 - - 



2 -^OOi 



8 -^ooi 



+ .... 



1,2) 



Using the Droplet model [22], neutron skin ti, reads as 



ti 



-Jo 



-1/3 



(7) 



(8) 



Here Tq is 1.14 fm. the value of nuclear symmetric energy coefficient J=32.65 MeV and 
Ci = 3e^/5ro = 0.757895MeV. The neutron skin stiffness coefficient Q was taken to be 
35.4MeV. 



The fragmentation potential Vjiir]) (at R — Ra — Ct{— Ci + C2))in (2) is calculated 
simply as the sum of the Coulomb interaction, the nuclear proximity potential [23]with 
new isospin dependent radii and the ground state binding energies of two nuclei. 

Via, ry) = - V B{A, Z,) + + Vp, (9) 



i=l 



Ct 



The proximity potential between two nuclei is defined as 

Vp = 47rC76$(0 



(10) 



here 7 is the nuclear surface tension coefficient, C determines the distance between two 
points of the surfaces, evaluated at the point of closest approach using eq.(4) and $(^) is 
the universal function, since it depends only on the distance between two nuclei, and is 
given as 

' -0.5(^-2.54)2-0.0852(^-2.54)3, 
for C < 1-2511 
-3.437ea;p(-C/0.75). 
far ^ > 1.2511 

4 



^(0 



(11) 



Here, ^= s/b, i.e s in units of b, with the separation distance s=R — Ci — C2. b is the 
diffuseness of the nuclear surface, given by 

h= [7r/2\/31n9] (12) 

J tlO-90 

where tio-go is the thickness of the surface in which the density profile changes from 90% 
to 10%. The 7 is the specific nuclear surface tension, given by 

2" 



7 = 0.9517 



1.7826 



N - Z 
A 



MeVfm~\ (13) 



In recent years many more microscopic potentials are available that takes care various 
aspects such as overestimation of fusion barrier in original proximity potential, isospin 
effects. A comparison is also available between all models [24]. The binding energy 
B{Ai, Zi) (i=l,2) of any nucleus of mass number A and atomic number Z was obtained 
from a phenomenological search and was given by a more genaralized BW formula [25] 
used and found good in agreement with experimental results by many others for drp- 
line to superheavy nuclei [26]. Thus, shell effects are also contained in our calculations 
in addition to the isospin effects for all the normal to neutron/proton rich nuclei. The 
momentum dependent potentials and symmtry energy potential which are found to have 
drastic effect at higher densities will not affect decay studies, since these happens at lower 
tale of the density [27,28]. Here in Eq. (4), the Coulomb and proximity potentials are 
for spherical nuclei, and charges Zi and Z2 in (4) are fixed by minimizing the potential in 
rjz coordinate. The mass parameters Bj^^{r)), representing the kinetic energy part in Eq. 
(2), are the classical hydrodynamical masses of Kroger and Scheid [29]. 

The WKB tunnelling probability, calculated is P = PiPj, with 

p. = exp^l [ ' {2ix[V{R) - V{Ri)]y/^dR\ (14) 

Pb = exp[-\ [ \2fx[V{R) - Q]y/^dR]. (15) 

'^J Ri 

These integrals are solved analytically [16] for the second turning point, defined by 
V{Rb) = Q- value for the ground-state decay. 

The assault frequency uq in (1) is given simply as 

= {2E2/i^y/yRo, (16) 
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Figure 1: (Color online) (a) logioTy^i^) for a-decay series of ^^"^117 nucleus. (b) logiQTy2{s) 
for a-decay series of '^^'^\YJ nucleus, calculated on the basis of the Isospin Cluster Model 
(ICM), compared with the experimental results, preformed cluster model (PCM) and ICM 
as a function of calculations, plotted as a function of parent mass number, in these 

figures. 



with E2 = {Ai/A)Q, the kinetic energy of the lighter fragment, for the Q-value shared 
between the two products as inverse of their masses. 



3 Calculations and results 

In Fig.l two a-decay chains of the Z—117 isotopes, ^^^117 and ^^^117 are shown for the 
alpha decay half-hves and Q-value calculations. A comparision of PCM and ICM cal- 
culations with experimental results are shown. The Isospin Cluster Model calculations 
are comparable and are in good agreement with experimental results. For both the al- 
pha decay chains ICM calculations also performed as a function of Q^^p. The results 
ICM{Q^^^) are resonably well for the Z — 107, 109 nuclei in comparision to ICM with 
the others in ^^^117 a-decay chain. Where Q^'^'^- value is calculated from the BW mass 
formula suggested by samanta et.al. [25]. Two a-decay chains are calculated within the 
ICM model framework for the four different Q- values i.e., Q^^^, Q^"^, QR^FiNLs)^ 
QiCM ^Yie effect of the Q-value on the alpha decay half-hves. Which can be seen 
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Figure 2: (Color online) (a) The Q^-energy and (b) the half-life times logx^T'^i^i^s) for 



a-decay series of ^^^117 nucleus, calculated on the basis of Isospin Cluster Model (ICM), 
compared with the experimental results [9], Macro- Microscopic [10], RMF(NL3 parame- 
ter) [11], RMF( FSU Gold forces parameter) [12] and ICM as a function of Qf^^'*- calcu- 
lations, plotted as a function of parent mass number, in these figures. 
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Figure 3: (Color online) Same as for Fig. 2, but for the ^^"^117 nucleus. 
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in the Fig. 2 and 3 clearly. Further the ICM Q-values Q are compared with the 
experimental (Q^^^), Q'^^^ , qR^^'C^ls) v^i^^gg fo]; both the alpha decay chains. In table 
1, Fig. 2 and 3 the ICM model calculations for the alpha decay half- lives are compared 
with the experimental [9], Macro-Microscopic calculations (M-M)[10], RMF model calcu- 
lations using NL3 [11] and FSU Gold forces [12] parameters. In Fig. 2b the calculated 
ICM half-lives results are comparable and reasonably well with experimental and other 
theoretical model calculations. The ICM half-lives as a function of Q^^p^- are also studied 
and shown in figure. 

Similarly, in Fig. 3, a comparision of different Q-values and half-lives are shown for the 
^^^117. The calculated ICM and ICM function of {Q^^^^')j half-lives are in good agreement 
with experimental results except to the ^^^111 [Z = lll.A^ = 171)nucleus, where a 
noticeable maximum half-life is found for the nucleus indicates the shell effects(N=172). 
The decay studies shows that half-hves of the alpha decay as well as cluster decay works 
as a tool in nuclear structure physics to show the persence of shell effects of the parents 
as well as of daughter nuclei. A higher value of the half-hfe indicates the presence of 
shell stabihzed parent nucleus, whereas a comparatively low value of half-hfe tells the 
same about the daughter and cluster nuclei [14,15]. The RMF(NL3) calculation shows 
the same trends and results except to ^^°115 {Z = 115, iV = 175) nucleus, with the ICM 
calculations. The other RMF calculation using the SFU Gold parameter also follows the 
same trends as the ICM(g^^P* ) results. 

4 Discussion and summary 

The Isospin Cluster Model (ICM) is used for the calculations of the two ^^^117 and 
^^^117 alpha-decay chains. The calculated alpha decay half-hves are compared with the 
experimental results and the other theoretical models, which are in good agreement with 
experimental results and comparable to the other theoretical model calculations (M-M, 
RMF). The a-decay calculations are performed for the two decay chains of 293,294-1^27 nuclei 
keeping the parents and daughter in spherical nuclear shape and transtion of alpha decays 
are considered from ground state of the parents to the ground state of daughter nuclei. A 
maximum half-life at -^^^-^lll {Z = 111,A^ = 171) in ICM{Q^^p^-) calculation is observed 
indicates the shell effects of the parent nucleus, near to the proposed shell closure (N=172) 
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by other RMF model calculations. 
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Table 1: Comparison of calculated Isospin Cluster Model(ICM) «-decay half-lives with 
experimental results, Macroscopic-Microscopic (M-M) [10] and with RMF model calcu- 
lations (using NL3 [11] and FSU Gold forces[12]parameters). Q^^^ is calculated using 
binding energies from BWMF [25]. 



Parent gf^P[9] Q^^[10] Q^^^[2h\ Exp.[9] MM[10] ICM ICM (gf^f- 





(MeV) 


(MeV) 


(MeV) 


LogioTii2{s) 


LogioTii2{s) 


LogioTii2{s) 


LogioTi/2' 


293117 


11.183 


11.42 


11.699 


-1.837 


-2.159 


-0.797 


0.059 


289115 


10.45 


10.63 


11.100 


-0.654 


-0.817 


-0.149 


1.041 


285113 


9.88 


10.10 


10.495 


0.738 


-0.079 


0.622 


1.87 


294117 


10.97 


11.15 


11.56 


-1.109 


-1.506 


-0.402 


0.269 


290115 


10.10 


10.37 


10.95 


1.797 


-0.159 


0.147 


1.817 


286113 


9.77 


9.70 


10.34 


1.293 


-1.955 


0.941 


2.139 


282111 


9.13 


9.57 


9.73 


-0.289 


0.749 


1.805 


3.202 


278109 


9.67 


9.27 


9.11 


0.882 


0.955 


2.821 


1.55 


274107 


8.93 


8.55 


8.47 


1.732 


2.488 


4.06 


2.89 
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